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Abstract 

The problem of comparison of the stationary axisymmetric vacuum 
solutions obtained within the framework of exact and approximate ap- 
proaches for the description of the same general relativistic systems is 
considered. We suggest two ways of carrying out such comparison: (i) 
through the calculation of the Ernst complex potential associated with 
the approximate solution whose form on the symmetry axis is subse- 
quently used for the identification of the exact solution possessing the 
same multipole structure, and (ii) the generation of approximate so- 
lutions from exact ones by expanding the latter in series of powers of 
a small parameter. The central result of our paper is the derivation 
of the correct approximate analogues of the double-Kerr solution pos- 
sessing the physically meaningful equilibrium configurations. We also 
show that the interpretation of an approximate solution originally at- 
tributed to it on the basis of some general physical suppositions may 
not coincide with its true nature established with the aid of a more 
accurate technique. 



PACS numbers: 0420 
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1. Introduction 

Approximation methods had long proved to be an efficient tool for treating 
specific general relativistic problems. For many years their use was predom- 
inant since only a restricted number of physically realistic exact solutions of 
Einstein's equations was available before the development of modern solu- 
tion generating techniques [1-7]. But even now, when many exact solutions 
have already been constructed, approximation methods are frequently used 
for attacking some interesting problems not yet solved exactly. Among the 
relatively recent successful applications of the approximation schemes one 
can mention the analysis of the frame-dragging effect in the field of a static 
charged massive magnetic dipole jS] with the aid of the Martin-Pritchett 
approximate solution jU] , or a study of the equilibrium problem of two static 
charged masses using a test-particle approach [TU] . The exact solutions which 
are in accord with the above approximate findings were constructed later on 

In spite of a rather peaceful coexistence of exact and approximation ap- 
proaches it appears that there does exist a problem upon which the two 
approaches had produced some not fully concurring results - the vacuum 
equilibrium problem of two spinning particles. There is only agreement of 
the exact results with the well-known paper of Wald [13 in which he proved 
within the framework of a test-particle approach that the spin-spin interac- 
tion cannot balance the gravitational attraction of two point-like particles. 
In his effort to develop a far-reaching successive approximation procedure 
which would go beyond the paper ^3], Bonnor derived an approximate so- 
lution which he considered applicable to any two spinning uncharged 
particles, black holes or superextreme objects. In that solution, however, the 
balance was not achievable too, this time producing an apparent contradic- 
tion with the known exact results. Indeed, it is well known that equilibrium 
states of two superextreme particles do exist in the double-Kerr solution first 
constructed by Kramer and Neugebauer The general analytic solution 
of the extended double-Kerr equilibrium problem was given in and 
in a simple relation was found between the coordinate distance at which 
the balance occurs and the masses and angular momenta of the balancing 
Kerr constituents. 

1 Equilibrium states of two identical superextreme Kerr constituents were first obtained 
and analyzed by Dietz and Hoenselaers |18| . 
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This gives rise to an important question: how can one compare exact and 
approximate solutions aimed at describing the same physical situations? It is 
very important to have an answer to this question for being able to interpret 
approximate solutions correctly, establish limits of their validity, and give a 
recipe for the rectification of an approximation scheme which for some rea- 
sons has turned out to be incompatible with an exact model. In this paper 
we shall develop an effective approach to the comparison of the exact and 
approximate solutions based on the following mutually complementary pro- 
cedures: (i) identification of the exact analogue of an approximate solution 
through the axis expression of the Ernst potential j20] of the approximate 
solution, and (ii) generation of the approximate solution from an exact one 
by expanding the Ernst potential of the exact solution in a series of pow- 
ers of a small parameter with a subsequent comparison of the result with 
the known approximate solution. Whereas the procedure (i) provides one 
with a precise interpretation of the approximate solution (which may not 
coincide with the interpretation initially attributed to it), the procedure (ii) 
allows the construction of the approximate solution best matching to the 
exact one. To illustrate our approach, we shall apply it to the classical prob- 
lem of two spinning particles, taking as reference points the recent Bonnor's 
approximate solution [T3j and the exact double-Kerr solution in its extended 
form ^7]. Note that the solution was claimed to be appropriate for the 
approximate description of two Kerr particles, but we will show that this 
affirmation is erroneous. 

The paper is organized as follows. In Section 2 we derive the Ernst 
potential of the Bonnor solution and give the interpretation of the latter 
by means of the corresponding exact solution. In Section 3 we construct 
the approximate double-Kerr solution in the subextreme case not admitting 
equilibrium configurations. Section 4 deals with the approximate double- 
Kerr solutions involving one or two superextreme particles for which the 
corresponding equilibrium problems can be solved in full analogy with the 
exact solution. In Section 5 we establish the relation of the parameters 
entering approximate solutions to the individual Komar masses and angular 
momenta of the particles. Section 6 presents a discussion of the results 
obtained. Concluding remarks are contained in Section 7. 
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2. Interpretation of Bonnor's approximate solution 



In this section we will show how a precise physical interpretation of an ap- 
proximate axisymmetric solution can be established via an appropriately 
identified exact solution whose properties can be ascribed to the approxi- 
mate solution. Though the existence of a solid basis for finding such an 
exact analogue may look questionable at first sight, this is fortunately not 
the case. The well-known Sibgatullin's integral method of the construction 
of exact solutions [7j permits one to generate axisymmetric spacetimes from 
the axis expressions of their Ernst complex potentials [20] which are used 
as the initial data in the generation process. This means that all one needs 
for having at hand an exact analogue of a known approximate solution is to 
establish the form of the Ernst potential on the symmetry axis of the lat- 
ter; Sibgatullin's method will then provide one with the corresponding exact 
solution possessing the same axis expression of the Ernst potential, and its 
properties can be studied with the aid of the well established procedures. 

This scheme can be used for establishing the physical interpretation of 
the approximate vacuum solution for two spinning particles ^3]. Since this 
solution was obtained in a somewhat special form, we first of all have to 
rewrite it using the canonical axisymmetric line element 2 

ds 2 = f-^e^idp 2 + dz 2 ) + p 2 dip 2 } - f(dt - ood V ) 2 , (1) 

where the metric coefficients /, 7 and u depend on the coordinates p, z only; 
<p and t are the athimuthal angle and time, respectively. 

Observing that the line element (1) of [Hj can be cast, after changing the 
signature, into the form 

ds 2 = f- l [ e »(dr 2 + dz 2 )+r 2 d6 2 }- f{dt-nf- l d9) 2 

- r\r 2 -n 2 -lf)de 2 , (2) 

we obtain formulae relating notations of the paper and ours: 

r = p, z = z, v = 27, n = fu, r 2 — n 2 — If = 0. (3) 

Now, bearing in mind (jSJ) and using formulae (14)-(20) of Ref. |14j . we 
are able to rewrite the Bonnor solution in its standard form: 

t - ^'(f +fM (f +f) 2 

2 Throughout the paper units are used in which Newton's gravitational constant G and 
the speed of light c are equal to unity. 
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In the above formulae we have left Bonnor's notations for the masses mi, m 2 
but instead of his hi, h 2 we have used ours ji = —hi, j 2 = —h 2 as angular 
momenta of the particles; the parameter b defines location of the particles 
on the symmetry axis (the points z = ±6). We have introduced an auxiliary 
nondimensional small parameter I by the formulae 



/ = /(°) + I fW + / 2 /( 2 ), v = V W + l^m + /V 2 ), 



n 



(5) 



(f( l \ are defined in index i denoting the order of approxima- 

tion), and it helps to control the order of approximation in Bonnor's scheme 
(one can always put I = 1 at the end). Note that the metric functions P|) 
cannot contain terms with l n , n > 2, because the solution is given in the sec- 
ond approximation. Mention also that the masses rrii and angular momenta 
ji are quantities of the same (first) order in I. 

The metric coefficients /, 7, to defined by (@J) are approximate solutions to 
Einstein's equations in the stationary axisymmetric case (see, e.g., Section 7.1 

of EH) 



/A/=(V/) 2 -p- 2 / 4 (V^ 
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V(p- 2 / 2 Vc) = o, 

7,p = V 2 [/j-/j-p- 2 /vi-^)], 



i, z = ^ P rvj, z -p- A r^x z ), (6) 

where a comma in subindices denotes partial differentiation, and the opera- 
tors A, V have the form 

d 2 1 d d 2 -* _ d _ d 
op 2 pop oz 2 op oz 

( P o and Zq are unit vectors). 

The first two equations in (0) can be rewritten as a single Ernst equation 
for a complex potential £ 



{£ + £)A£ = 2(V£) 2 , (8) 

a bar over a symbol denoting complex conjugation. 

The relation of £ to the functions / and uj is the following: 

e = f + in, n, p = p-W*. n^-p-'/Vp. (9) 

Sibgatullin's method [7] provides a straightforward procedure for the con- 
struction of exact solutions of equation (JHJ) corresponding to a given function 
e(^) = £(p = 0,z) representing the axis value of the Ernst potential £. 
Therefore, in order to identify the exact analogue of Bonnor's approximate 
solution (JU), we have to find from (UJ) and Q the axis value es{z) of the cor- 
responding Ernst potential in the second approximation. Having obtained 
£b{z), we shall be able to point out what exact solution corresponds to this 
function and, besides, to use cb{z) for the calculation of the multipole mo- 
ments of both the approximate and exact solutions with the aid of the Fodor 
et al. procedure [22] • 

The integration of equations © for Q, with / and to defined in 03} , can 
be readily performed using the Mathematica computer programme [211] • The 
resulting expression, multiplied by i, should then be added to /, yielding the 
desired Ernst potential of the Bonnor solution 

£ B = l- 21 M + 2l 2 (M 2 - R 2 ) - 2ilR 

2 / z(mij 2 + m 2 ji)(p 2 + z 2 - b 2 ) + b(m 1 j 2 - m 2 ji)(p 2 - z 2 + b 2 ) 

I mm 



. mdi(z-b) m 2 j2(z + by 
H r 



K x K 2 J 

M =ik + 7i b =t +j t < 10 » 

On the upper part of the symmetry axis, p = 0, z > b, the potential 
obtained takes the form 

„ , v 2/iV 2l 2 N 2 

e B (z) = £ B {p = 0,z) = 1 - —. — + 



y z 2 -b 2 ) 2 (z 2 -b 2 ) 4 ' 
N = (z 2 — b 2 )[z(rrii + m 2 ) + b(m 1 — m 2 )] 

+ i(ji+j2)(z 2 + b 2 ) + 2ibz(j 1 -j 2 ). (11) 

Expression (|TT| is the ratio of two polynomials of the eighth order in z, 
so the exact solution corresponding to this axis data necessarily falls into the 
family of 2iV-soliton solutions discussed in detail in [21]. Moreover, as the 
denominator of es{z) is the real function of z, it should be a special member 
of the N = 4 subfamily of the soliton solutions which in the general case 
involves 16 real parameters, contrary to only 5 parameters of the axis data 

dm. 

To see what physical situation describes the exact solution whose Ernst 
potential on the symmetry axis has the form one has to recall that 

after the application of Sibgatullin's method to the data (fTTjl the resulting 



solution will contain the functions = y p 2 + (z — «j) 2 , where are roots 
of the algebraic equation 

e B {z) + e B {z) = 0. (12) 

Apparently, this equation has 8 roots which can assume real values or occur 
in complex conjugate pairs, hence the solution will describe special configu- 
rations of 4 spinning objects because each pair of aj determines one sub- or 
one superextreme constituent. 

For a more precise analysis we need the explicit form of Obviously, 
equation (J12J) cannot be solved exactly, but its approximate roots up to the 
second order in I are available. These should be searched for in the form 
±b + Z\l x l 2 + z 2 l + z^l 2 + zj 2 , the constant coefficients zi, z 2 , z 3 , z 4 to be 
found from (|12|) . Then one gets the following 8 roots of equation (|12|) : 



Oil 



+ Jnl 1 ' 2 + A^' 2 + A 12 l 2 
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a 2 = b- ijnl l/2 + lA^' 2 + A 12 l 2 , 
a 3 = b + iJJll 112 - iA^ 2 + A 12 l 2 , 
a 4 = b - ^l 1/2 - A x l zt2 + A 12 l 2 , 
a 5 = -b + J^l 1 ' 2 + A 2 l 3 / 2 - A l2 l 2 , 
a 6 = -b- i^j~ 2 l 1/2 + iA 2 l 3 / 2 - A 12 l 2 , 

a 7 = -b + i^]~ 2 l 1/2 -iA 2 l 3/2 - A 12 l 2 , (13) 

« 8 = -b - VS 1/2 - AJ 3 ' 2 - A 12 l 2 } 

_ 2b 2 m\ + jij 2 _ 2b 2 m\ + j 1 j 2 2b 2 m 1 m 2 - j x j 2 

86 2 v^ ' 2 86 2 v^ ' 86 3 

which means that independently of the sign of j\ and j 2 these as determine 
a system formed by two sub- and two superextreme constituents. More pre- 
cisely, we have a system of two separated compound objects, each of which is 
composed of the overlapping sub- and superextreme constituents. One easily 
arrives at this conclusion by analyzing for instance the expressions for a%, 
a 2 , «3, «4, of which a± and define a subextreme constituent (assuming 
ji > 0), whereas a 2 and a 3 define a superextreme constituent; these two 
overlapping constituents form a compound spinning object. Another com- 
pound object is formed by a sub- and a superextreme constituents defined 
by the constants a 5 , a e , a 7 , a 8 (see Fig. 1). 

Therefore, the physical interpretation of the Bonnor solution (0J) supplied 
by the corresponding exact solution is the following: it represents the second 
approximation to a special case of the N = 4 subfamily of multi-soliton exact 
solutions describing two compound objects located on the symmetry axis and 
kept apart by a strut in between. Such systems of two compound objects have 
been considered in the paper 

It should be also pointed out that Bonnor's solution cannot be consid- 
ered as an approximation to the well-known double-Kerr spacetime [To] . 
Although the general iV = 4 soliton metric (the quadruple-Kerr solution) 
possesses 16 arbitrary real parameters and permits in principle the reduction 
to the simpler N = 2 double-Kerr case, such a reduction is not possible for 
as defined in (|T3*|) because of the insufficient number of constants required 
for performing the desired reduction. Correct approximations to the double- 
Kerr solution will be obtained in the next two sections. 

To conclude this part, below we write down the first five complex mul- 
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tipole moments of the Bonnor solution calculated from (jlljl with the aid of 
the Fodor et al. procedure 

P = l( mi +m 2 ), Pi = l[b(m 1 -m 2 ) + +J2)], 
P 2 = U[b(mi + m 2 ) + 2i{j 1 - j 2 )], 
P 3 = 6 2 /[6(m! - m 2 ) + 3i(j'i +J2)], 

P 4 = 6 3 Z[6(772! + m 2 ) + 4i(j! - j 2 )]. (14) 

The form of these multipoles confirms the interpretation of m\, m 2 as 
masses, and ji, j 2 as angular momenta. However, a more subtle analysis is 
needed to see whether m; and ji coincide with the individual Komar masses 
and angular momenta [26J of the particles (see our Section 5). 

3. Approximate double— Kerr solution in the black— hole case 

When a problem has an exact solution of a complicated form, one may think 
about the conversion of that solution into an approximate one under some 
simplifying assumptions. The generation of approximate solutions from the 
already known exact solutions has obvious advantages compared to the de- 
velopment of approximation schemes exclusively on the basis of some general 
suppositions and physical intuition, first of all because the approximation in 
the former case will not contradict the existing exact results but will only 
permit their presentation in a simpler form. On the other hand, there is 
always a danger that an approximate method not orientated at the exact 
solution will lead to a result contradicting the exact findings and, besides, 
will be unable to identify the origin of the contradiction. 

Usually the exact solutions constructed with the aid of the modern solu- 
tion generating techniques, especially those describing the many-body con- 
figurations, involve the parameter sets most adjusted to the mathematical 
structure of the solutions but which may look "unphysical" to non-experts 
of the generation methods. At the same time, a rather sophisticated relation 
of the parameters of the generated solutions to the physical characteristics 
of spacetimes they describe certainly may course difficulties for introducing 
the simplifying physical assumptions for a possible approximate analysis of 
the generated solutions. 

The aim of this section and of the next one is to demonstrate that within 
the framework of Sibgatullin's method it is possible not only to construct ex- 
act solutions but also develop a consistent approach to the analysis of their 
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physical properties on the basis of reliable approximation schemes. Below we 
shall obtain approximate analogues of the well-known double-Kerr solution 
15J making use of the unique opportunities offered by Sibgatullin's method 
- the possibility of the unified treatment of the sub- and superextreme cases 
of spinning particles and a very clear relations of all the parameters of the 
solution to the axis data and, consequently, to the relativistic multipole mo- 
ments. The Ernst potential of the extended double-Kerr solution has the 
form J7| 



£ — a _ p) A — ^ ^ij r i r j> F — ^ h ' ' 
A 1 l<i<i<4 i=l 

^ij = (- l T +J ( a i ~ <Xj)(an> - aji)XiXj, (i',f ^ i < f), 
v% = {-l)\a v - ajt){a v - a k >)(a r - a k ,)X h 
(*',/, i' < f < k') 

(a* - /3i)(«i - (3 2 ) 



Xi = Tr, — wxu — aT' r i = yp 2 + ( z -^) 2 , (15) 

[ati - pi){ai - P2) 

The parameters entering (JT5j) are (3\, @2 which can assume arbitrary complex 
values, and aj, i = 1,2, 3, 4, which can take arbitrary real values or occur in 
complex conjugate pairs (without loss of generality we can assign to aj the 
following order: Reai > Re«2 > Re«3 > Rec^). On the symmetry axis £ is 
defined by the expression 

«(«) = i + A + ^&- (16) 

where e\ and e 2 are two arbitrary complex constants. 

Note that are roots of Sibgatullin's algebraic condition 

e(z) + e(z) = 0, (17) 

which means that the following relation holds: 

2+ePv + ^u= 2 2nti( :~ a " ) , v as) 



i=i 



whence one obtains the connection between the parameters e; and a n 

= 2nLi(A-«») = 2nUM-<*n) (m 

1 ah - ihM - Pirn - fay 2 - Pi)(jh - Pi)(P2 - fay ' 
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In the papers [T71 HH] we presented the general solution to the equilib- 
rium problem of two Kerr particles and found a simple relation between the 
coordinate distance at which equilibrium takes place and individual Komar 
masses and angular momenta of the balancing constituents. The question 
which now would be of interest to answer in view of the paper ^3] is whether 
the double-Kerr equilibrium problem can be solved with the aid of some 
approximation procedure? Below we are going to show that the answer to 
this question is 'yes' and it is supported by an approximation scheme which 
is able to reproduce the exact results. 

In order to develop the desired physically meaningful approximation pro- 
cedure for the double-Kerr solution it is likely first of all to rewrite the formu- 
lae (|15|) in a more "physical" parameter set which would be mathematically 
equivalent to the parameter set used in (|15j). The reparametrization of the 
axis data (fTBj) can be done in analogy with a single Kerr solution |2Zj whose 
axis expression of the potential £ = (1 — £)/(! + £) has the form 3 

. . 1 — e(z) m . . 

1 + e{z) z — la 

where m is the total mass and a is the angular momentum per unit mass of 
the Kerr source. 

In the case of the double-Kerr solution we simply generalize (|2U|) to the 
expression 

m l + iv m 2 -iv 
x \ z ) = ; : 1 ; : — > (21) 

Z — — Idi Z + — ICL2 

taking mi, rri2 as the "masses" of the Kerr particles, and ai, a 2 as their 
"angular momenta per unit mass" ; the parameter b denotes displacements of 
the particles from the origin of coordinates along the z-axis, and v represents 
the angular momentum of the part of the symmetry axis separating the 
particles. 

From (j21j) we obtain the corresponding expression of the Ernst potential 

£ 

. . 1 — x(z) e_ 

e(z) ^ 



1 + x(z) e+ 



3 The potential £ = {£ — !)/{£ + 1) was first introduced by Ernst in his famous paper 
plj . Changing the sign in the original expression is explained by the successful use of 
the modified potential in the Fodor et al. procedure [221 for the calculation of multipole 
moments. 
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e± = (z — b — ia\){z + b — ia 2 ) 

± {mi + iv)(z + b — ia 2 ) ± im 2 — iu)(z — b — iai). (22) 

From (}2*T]) the expressions of the total mass M and total angular momen- 
tum J of the system can be readily found with the aid of the Fodor et al. 
procedure: 

M = mi + m 2 , J = micii + m 2 a 2 + Ibv. (23) 

This supports our physical interpretation of the parameters, even though 
their precise relation to the individual Komar masses and angular momenta 
will be established later on in Section 5. It is worth pointing out that the 
parametrization (J22j) of the Ernst potential involves 6 arbitrary real parame- 
ters and it has one advantage compared to the axis expression which is 
the absence of the unphysical NUT parameter that guarantees the asymptotic 
flatness of the solution from the very beginning. Up to this NUT parame- 
ter, the two parametrizations are totally equivalent in a sense that both of 
them represent all possible stationary asymptotically flat configurations of 
two Kerr particles. 

Now it is possible to develop an approximation scheme as follows. We can 
assume the values of the parameters mi, m 2 , ax, a 2 , v to be small compared 
to the values of b, which is equivalent to considering them proportional to 
some auxiliary non-dimensional small parameter I. Then an approximate 
expression of the potential £ can be found as an expansion in powers of / 
up to the desired order. Technically this consists in (a) finding approximate 
expressions for ot n from equation (fTTj). {b) identifying the expansions for 
X n from the formulae (j!5|) . (}2*2*|) in which the approximate expressions for 
ot n should be used, (c) obtaining the expansions for r n corresponding to 
approximate as, and finally (d) the expansion of the potential £ in powers of 
/ after the substitution of a n , X n , r n previously obtained into the formulae 
(JT3J). Once the potential £ is found, the approximate expressions for the 
metric coefficients uj and 7 can be derived from equations Q and © (note 
that in spite of the availability of the exact expressions for uj and 7 ^E] it 
looks more attractive to carry out additional integration instead of working 
out the expansions of the known formulae). 

We have to distinguish between the following three major approximation 
schemes which must be treated separately. 

Scheme 1. All five parameters m\, m 2 , a 1; a 2 , v are assumed to be small 
compared to b. In this case, making the substitutions m ; — > Imi, a { — > la i: 



12 



v — > Iv in (|22j) . we see that Sibgatullin's condition (|17|) at the zeroth order 
in I yields the equation 

( z 2 - b 2 ) 2 = 0, (24) 

whose double roots +b and —b should be taken as the zeroth approximation 
for the parameters c^: a>i = a 2 = b, a 3 = a 4 = —b (see Fig. 2a). 

Scheme 2. In this scheme, mi, m 2 , v are small parameters compared to 
b, while ai and a 2 are of the zeroth order in I, i.e., 'comparable' with b. After 
the substitutions m^ — > Irrii, v — > lv, the zeroth order in I of equation (fT7j) 
takes the form 

(z 2 - 2bz + b 2 + al)(z 2 + 2bz + b 2 + a 2 2 ) = 0, (25) 

and the roots of this equation provide us with the zeroth approximation of 
af «i = b — iai, «2 = b + idi, a 3 = —b — ia 2 , a 4 = —b + ia 2 (Fig. 2b). 

Scheme 3. In this scheme we assume that only one of the constants a\, a 2 
is comparable with b, say, a 2 , whereas m 1; m 2 , a 1; z/ are small parameters. 
Equation (|T7|) in this case leads to 

(z - b) 2 (z 2 + 26^ + 6 2 + a 2 ) = 0, (26) 

giving the following zero-order approximations for c^: a>i = a 2 = 6, «3 = 
—6 — m 2 , a 4 = — 6 + ia 2 (Fig. 2c). 

Remarkably, the schemes 1-3 correspond to the configurations composed 
of two subextreme constituents, of two superextreme constituents, and of 
one sub- and one superextreme constituents of the extended exact double- 
Kerr solution, respectively (we remind that in the exact solution a pair 
of real-valued as defines a subextreme constituent, and a pair of complex 
conjugate as defines a superextreme Kerr constituent). Mention also that 
only the first scheme can be considered as a sort of a "point-like-particles 
approximation" since, as will be seen later on, it involves two real dis- 
tances R\ = \J p 2 + (z — b) 2 and R 2 = J p 2 + (z + b) 2 , whereas the remaining 
schemes involve, respectively, four and three different square roots Ri due to 
the presence of two (scheme 2) and one (scheme 3) superextreme objects 
whose absolute values of the angular momenta per unit mass exceed the 
individual masses. 

Below we shall work out in detail the approximation scheme 1; schemes 2 
and 3 will be considered in the next section. 

Approximate double-Kerr solution via scheme 1. Without loss of general- 
ity we can set b = 1 in the previous formulae (this simply means the rescaling 
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of the parameters and coordinates) for simplifying a bit the expressions we 
are going to obtain. Like the paper [Tlj . we shall be interested in the ex- 
pansion of £ and corresponding metric functions in powers of I up to I 2 (the 
second approximation). Then we have to calculate «j up to I 3 for obtaining 
the correct expressions for Xi in the second order of I. After introducing the 
factor I into the parameters rrii, cii, u, equation (J 17)) takes the form 

(z 2 -l) 2 - l 2 [(m\-a\)(z + l) 2 + (m 2 2 - a 2 2 )(z - l) 2 

+ 2m 1 m 2 (z 2 - 1) + 4z/ 2 ] + 2vl z (m l + m 2 ) [ai{z + 1) 

— a 2 (z — 1)] — / 4 [m 2 a 2 + ^2 a i + ^ 2 ( a i — a 2) 2 

+ aia 2 (2m 1 m 2 — aia 2 )] = 0. (27) 

The roots of this equation can be searched for in the form ±1 + z\l + 
Z2I 2 + z 3 l 3 , the coefficients Zi to be determined from (|27|). The resulting 
expressions for ctj are the following (in view of the cumbersome coefficients 
z 3 we write out as only up to the terms I 2 , but it should be remembered that 
the coefficients z 3 cannot be omitted during the calculations): 

I 2 I 2 a\v[m\ + 777.2) 



a.\ = 1 + I J m 2 + v 2 — a 2 H m^mi — v — 



2 y Wm 2 + z/ 2 — a 2 



Z 2 / 2 a\v[m\ + m 2 ) 



a 2 = 1 — / u m 2 + v 2 — a 2 H mim 2 — z/ + 



2 * - 'm 2 + v 2 — a 2 / 

I 2 I 2 CL2v{mi + 777,2) 



03 = — 1 + lym 2 , + z/ 2 — a 2 m 1 m 2 — z/ + 



m 2 + z^ 2 — a 2 



Z 2 / 2 a 2 z/(mx + m 2 ) 



a 4 = — 1 — Zym 2 + v 2 — a 2 mim 2 — v . (28) 



m 2 + v A — a 2 . 



The expansion of the constant objects X, from (JT5j) associated with the 
parameters can be most conveniently found with the aid of the formula 



(29) 



where e_i_ is defined in 
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The expansion of Xi should be performed up to the second power of Z; 
this yields (to avoid the complicated expressions, we give explicitly only the 
terms linear in Z) 



imf + v l — af + iai I ( 

Xi = — — — - 2a x a 2 - iv{mi + m 2 ) 

mi + iv l[vri\ + iv) V 



—2%a2\jm{ + v l — a{ + 



mf + v 2 — a 2 , 



\lm\ + v 2 — a 2 — iai I ( 

X 2 = — — ; — T 2«ia 2 - iv{m\ + m 2 ) 

mi + 2(mi + ) y 

aiv(m 1 + m 2 ) \ 



+2ia 2 v m i +v — cl{ — 



m\ + v 2 — a 2 , 



'm 2 + z/ z — a 2 + za 2 / / 

X 3 = ; h — — 2aia 2 + «V(mi + m 2 ) 

m 2 — 21/ 2(m 2 — iv) y 

o- — 2 2 a 2 K m i+ m 2)\ 

— 2za 1 A/m 2 + i/ 2 — a 2 — 



m 2 + z/ 2 — a 2 , 



v' m 2 + ^ 2 -a2-^2 Z / . . . 

A 4 = ; h — — 2aia 2 + iv[m\ + m 2 J 

m 2 — 2(m 2 — zi/J y 

a 2 v(mi + m 2 ) 



+2iai y/mi + v 2 - a| + v = . (30) 

ym^ + v 2 - a\) 

Before passing to £ there still remains to find the expansion of = 



p 2 + (z — a;^) 2 . The calculation of with the aid of (|2*Hj) does not exhibit 
difficulty and leads to the expressions 



(z- 1)VX + ^ 2 - af /(m 2 + v 2 - a\){z - if 
n ~ Rl ~ l R, 1 [ 2^ 

z{m,\m2 — v 2 ) — m\ — mim 2 + a\ a\v[m\ + m 2 )(z — 1) 



+ 



2R i IRJml + v 2 -a 2 



t (z-l)y/ml + ^-al ({m 2 + v 2 - a\){ Z - if 

r 2 — -Kl + I — I 



Ri 1 2R\ 
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z[m\m 2 — v 2 ) — m\ — m\m 2 + a\ a 1 z/(m 1 + m 2 )(z — 1) 



2ft 



r 3 



R2-I- 



z + l)\/m 2 + v 2 — a 2 

ft 



2R iy /ml + v 2 -a\ 

[m\ + v 2 - a 2 2 )(z + l) 2 
2ft 



z(mim 2 — v 2 ) + m\ + mim 2 — a 2 , a 2 v{m\ + m 2 )(z + 1) 



2ft 



r 4 = R 2 + l 



> + l) A /m| + z/ 2 -a| 

ft 1 



2ft ym 2 , + v 2 -a 2 2 

'(ml + v 2 - a 2 2 )(z+ l) 2 
2ft 



z{m\m 2 — v 2 ) + m 2 + m\m 2 — a 2 , a 2 ^(mi + m 2 ){z + 1) 



2ft 



2R 2 Jm 2 , + v 2 -a 2 . 



,(31) 



where we have introduced 



ft = VP 2 + (z - I) 2 , ft = VP 2 + ( z + !) 2 - 



(32) 



We have realized all the necessary steps for eventually being able to ob- 
tain the form of £ in our approximation scheme. The substitution of the 
approximate values of Xj and r« into (fTK|) yields, after expanding £ in 
powers of /, 



£ 



1 - 21 



f mi + iv m 2 — iv 



V ft 



ft 



2l u 



m\ + iv m 2 — iv 



ft 



ft 



mi (mi + ii/) (z — 1) ia 2 [m 2 — iv)[z + 1) 



ft 



(33) 



For completeness we also give the expanded version of the potential £ 
associates with £: 



1-8 
l+£ 



I 



+ il 



mi + iv m 2 — iv 



ft 



ft 



2 , ai(mi + iv)(z — i) a 2 (m 2 — iv)(z + 1) 



■ (34) 



The knowledge of the potential £ is sufficient for obtaining the corre- 
sponding metric functions /, w, 7 in p. The function / is simply the real 
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part of £, hence it has the form 



/ 



+ 



1 - 21 



a\v{z 



■1) a 2 v(z+l) 



mi 

— + } , 

Hi R 2 J 



m 2 y 2 
— I — v 



1 



1 - 
R 2 . 



Rf 



Rl 



(35) 



The metric coefficients uj and 7 can be found by integrating equations (JSJ) 
and (jSJl, respectively. The function Q entering equations © is defined by the 
imaginary part of £. With this, the integration of gives 



UJ 



2lv 



Ri 



z + l 
R 2 



2V 



v{nii + m 2 ) 



2 / miai mia 2 \ 



i^(mi + m 2 )(2: 2 + p 2 - 1) 
-R1-R2 



(36) 



The substitution of (|33j). (j3^j) into equations © and the integration of 
the latter system yields the form of the remaining function 7: 



7 



2 



(m 1 m 2 - u 2 )(z 2 + p 2 - 1) 2 

— — niim 2 + v 

KiK 2 



/ mi + v l m% + z/ z 
P 2 I ^ + " 



(37) 



Mention that the integration constants in the above expressions for uj and 
7 are chosen in such a way that both functions vanish at spatial infinity. 

The behavior of uj and 7 on the symmetry axis is crucial for establishing 
whether equilibrium of two Kerr particles can be achieved due to balance of 
the gravitational attraction and spin-spin repulsion forces. We recall that 
for two particles to be in equilibrium, the functions uj and 7 should vanish 
on the parts of the z-axis outside and between the particles (the particles 
are situated at the points z = ±1). A simple inspection of the formulae ()36|). 
(I37|) shows that uj = 7 = when \z\ > 1. By demanding uj = 7 = on the 
part \z\ < 1 of the symmetry axis, we arrive at two balance conditions 



- 4vl[l + l(nii + m 2 ) 
(condition of vanishing uj), and 

— l 2 {niim 2 — v 2 ) = 







(38) 
(39) 
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(condition of vanishing 7). 

Bearing in mind that the parameter / is an auxiliary parameter which 
can be always put equal to 1 in the final formulae, we see that the system 
(J58Jl - (j3l?)l has no physically acceptable solutions since u — satisfying (J3~H|) 
implies vanishing of one of the masses mi, m.2. 

This result about the absence of equilibrium states in the approximate 
double-Kerr solution is consistent with our theorem ^7] on the absence of 
balance of two Kerr black holes with positive masses in the exact double-Kerr 
solution. 

We turn now to two other approximation schemes within which equilib- 
rium configurations of two spinning particles are available. 

4. Two approximation schemes leading to the balance of Kerr 
particles 

The approximation scheme 2 involves two superextreme Kerr particles. In 
the exact double-Kerr solution a superextreme constituent is defined by a cut 
joining two complex conjugate as. As we have seen in the previous section, 
by supposing the angular momenta per unit mass to be greater than the 
masses, already in the zeroth approximation we obtained two pairs of complex 
conjugate as as roots of equation (j2*3j) . Hence, in the approximation scheme 
too a superextreme particle is represented by a cut, and not by a point on 
the symmetry axis. This obviously will lead to the appearance in scheme 2 of 
four different square roots Ri, i = 1,4, involving complex quantities, unlike 
in scheme 1 where we had only two real distances R\ and R2 defining the 
location of the particles on the symmetry axis. 

Here we shall consider a slightly modified scheme 2 in which the parameter 
v will be assigned the order 2 in I (the case when v is proportional to /, not I 2 , 
does not lead to equilibrium configurations). The roots of equation (fTTj) after 
the substitution mi — > Imi, ni2 — > lni2, v — > l 2 v into f!22l) have been found 
to have the following form (technically they are obtainable as in scheme 1): 

/ m? mim 2 \ ,,y(mi+m2) 
a x = oL 2 = 1 - zai +il 2 [zr L + - + I 3 



2ai ai + a,2 + 2i J a,\ + a>2 + 2i 

a, = a, = -1 - ia 2 + { ^ + - l ^ m 2 + m i - (40) 

\2a 2 ai + a 2 — 2t J a\ + 02 — 1% 

We have found it more advantageous for schemes 2 and 3 to calculate the 
quantities IXt instead of X h that simplifies the subsequent expansion of £. 
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The resulting expressions for are (IX 2 and IX4 are written without their 
quadratic terms in / because of their complexity) 



IX X = il 



2 7711(0! - a 2 + 2%) 
2ai(ai + a 2 + 2z)' 



^ 2iai(ai + a 2 - 2i) ^la x v{a x + a 2 - 2i) 

2 mi(ai — a 2 — 2i) m\ {a\ — a 2 — 2i) ' 
2 m 2 (ai - a 2 + 2i) 

lJi-3 — — ll r, 

2a 2 (ai + a 2 — 2%) 
2ia 2 (ai + a 2 + 2i) 2a 2 is(a 1 + a 2 + 2i) 

IJi^ — -. r r I 7T. r— . V*l) 

m 2 \a\ — a 2 — 2i) m 2 {ai — a 2 — 2z) 
The expansion of leads to the following result: 

,z - 1+ia-i ( m? mim 2 \ 

n = i?x - i/ 2 — - -J- + - 1 



i?x \2ai a! + a 2 + 2z ; 
r 2 = R 2 + il 2 + 



i? 2 \2<2i ai + a 2 — 2^ 

. j2 2; + 1 + ia 2 / m 2 mim 2 
^3 = #3 - «t = o 1- 



R 3 \2a 2 a% + a 2 — 2% J 

„ , 9 z + 1 — ia 2 / m 2 , mim 2 \ , . 

r 4 = i? 4 + z/ 2 -M + — 4- , 42 

i?4 V 2a 2 ai + a 2 + 2i I 



where 



R 1 = R 2 = v /p 2 + (z-l + ia 1 ) 2 , R 3 = R 4 = ^p 2 + (z + 1 + m 2 ) 2 . (43) 
From (JT5|) and (}4T)j) -()43 j) we find the form of the Ernst potential £: 



mi m 2 \ 2 / 1 1 

1 I — iv I 

-R 2 -R4 / V R2 Ra 



(44) 



The corresponding potential £ is defined by a particularly simple expres- 
sion 

m l , m 2\ ■ ,2 ( 1 1 



^ + _ ^ (45) 

V i? 2 -R4 / V R 2 R4 J 

Remarkably, the formulae (jUj) and (j43j) do not contain the functions 
Ri and R 3 . However, all R iy i — 1,4, enter the expressions for the metric 
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functions /, uj, 7 (the procedure of the calculation of these functions is fully 
analogous to the one used for scheme 1) 



/ 



7 



l-l 



mi rn-i wi2 rri2 
— - H H H 

Ri R2 -R3 Ra 



2mi7«2 




+ z 2 



+ v 



ml m% 2mim 2 

— + — + — ^T 7 

a x 02 ai + <2,2 — 2« 

2; — 1 + iai z + 1 + ia 2 



+ 



2mim 2 



ai + a 2 + 2? 

m 2 (,2 2 + p 2 -2,2 + a 2 + i) 



#1 



m 2 ,^ 2 + p 2 + 2z 



R3 
+ 1) 



+ c.c. 



a 2 R 3 R 4 



— 2mim 2 



x 



= Z 2 



+ p 2 + aia 2 — 1 — iz{a\ — a 2 ) — i(ai + a 2 ) 



2 2 
mi mi 

— - H + 

4a? 4al 



(ai + a 2 - 2i)R 2 R 3 
mim 2 



<2iRiR 2 



+ c.c. 



+ 



mim 2 



(ai + a 2 - 2i) 2 (ai + a 2 + 2i) 2 



m 2 [(^ - l) 2 + p 2 + a 2 ] m 2 [(^ + l) 2 + p 2 + 



4ajRiR 2 



x 



V + p 2 + aia2 



Aa 2 2 R 3 R A 
1 + iz(ai - a 2 ) + i(ai + a 2 ) 



(ai + a 2 + 2i) 2 RiR A 



— mim 2 



+ c.c. 



(46) 



('c.c' denotes complex conjugation of the terms in brackets). Interestingly, 
the constant v does not enter the formula for 7. 

Let us analyze now the possibility of balance of the Kerr particles in this 
approximation scheme. By construction, the functions uj and 7 vanish on 
the part \z\ > 1 of the symmetry axis. Hence, to achieve the equilibrium of 
the particles, we only have to fulfil uj = and 7 = on the part \z\ < 1 
of the symmetry axis separating the particles. The condition uj = with 
Ri = R 2 = 1 — z — iai, R 3 = R A = 1 + z + ia 2 yields 



4Z 2 



' 2mim 2 {ai + a 2 ) 
v ( ai + a 2 ) 2 + 4 



0, 



(47) 
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and the condition 7 = leads to 



4m 1 m 2 l 2 [(a 1 + a 2 ) 2 - 4] 



0. 



(48) 



[( ai + a 2 )2 + 4]2 



The solution of (|4Tj) is 



v = 



2mim 2 (ai + a 2 ) 



(49) 



(ai + a 2 ) 2 + 4 



and it fixes the choice of v in this approximation scheme. Equation (|4*K|) is 
fulfilled when ai+a 2 = ±2. Recalling that all the parameters are "measured" 
in the units of b, we obtain after introducing b explicitly: 



and this is the relation between the distance at which equilibrium takes place 
and the angular momenta per unit mass of the particles. 

Scheme 3. Let us consider now the remaining approximation scheme. We 
will show that it also permits equilibrium configurations, this time between a 
sub- and a superextreme Kerr particles, as it occurs in the exact double-Kerr 
solution. 

The substitutions to be made in the axis data ()22|) are the following: 
nil Irrix, m 2 — > lm 2 , ai — > la%, v — > I 2 v. This means that the particle 
with subindex 1 represents a subextreme constituent, while the particle with 
subindex 2 is a superextreme constituent (a 2 > m 2 ). Mention that choosing 
v proportional to I 2 is a manner to introduce the inequality v 2 < a\. As 
before, we shall develop this scheme in the second approximation. 

Under the above assumptions equation (fTTjl has the following four roots 
(the coefficients at Z 3 are not given because of their complexity) 



2b 



ai + a 2 |, 



(50) 




a 4 = — 1 — ia 2 + I 



2 m 2 [2m 2 + m 2 (2mi + m 2 )] 



(51) 



2a 2 (a 2 — 2i) 
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Then, after expanding IX^ up to I 2 , we get 
IX X = 



(a 2 _ 2i)(m 1 + ia^ + Jm\ — a 2 ) 



(a 2 + 2i)(mi - ia x + \jm\ - a?) mf(a 2 + 2i) 2 



x 



imia 2 (2m 2 — 2a\ + mxm 2 ) 1 + 



+ zz/(a 2 + 4)(Jm\ — a 2 + ia x ) 



IX 2 
IX 3 

IX A 



X 1 {Jm\ 



m\ 



2 im 2 



2ia 2 2a 2 l ( im 2 a\ 



m 2 



+ 



ml 



— v 



2a 2 (m\ + 2a\ 



a 2 + 2i 
) Mm 2 a\a 2 v 



if 



m 2 2a 2 v 
2a>? m 3 



m 2 (a 2 + 2i) 2 m 2 i (a 2 + 2i) / 
The expansion of the functions r« yields 

I 2 



(z-l)y/ml-al 

tlx — I — h 



Ri 



2m\m 2 {z — 1) 
a? + 4 



2R X 
a x a 2 
m? — a? 



1 - 



'2 

r4 



i?4 — / 



m 2 (z + 1 + w 2 )[2m 2 + ia 2 (2m 1 + m 2 )] 
2a 2 (a 2 - 2i)i? 3 
2 m 2 (z + 1 — ia 2 )[2m 2 — ia 2 (2m\ + m 2 )] 



2a 2 (a 2 + 2i)i?4 



where we have introduced 
R x = Jp 2 + (z- 



l) 5 



— Ra 



p 2 + (z + 1 + ia 2 ) 2 . 



(52) 



(53) 



(54) 



The potentials £ and £ resulting from the formulae (jBT|) - (f5^|) and 
have the form 



£ 



2H^ + ^) + 2l 2 



mi m 2 \'- 
Ri R4 J 
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— IV 



1 1 



. -Ki -K4 



%m\a\{z — 1) 
miai(z — 1) 



R\ 



V 







(ri 


Ri) 



(55) 



The set of the metric functions /, u, 7 defined by ()55|) is given by the 
expressions 



./ 



I 



7? R R ' 

iXi 1X3 1X4 



+ 



mi 



2/ 



'm 2 (-2 + 1 + z'a 2 ) 
R~i 



1 



mi m 2 \ 2 
i?i i? 3 / 



m 2 [4m 2 + 03(47711 + m 2 )] 



+ 1/ 



a 2 (a 2 + 4) 

2(z — 1) 2 + 1 + ia 2 ,2 + 1 — ia 2 \ 2rriiaip 2 



Ri 



R* 



Ri 



RI 



m 2 [{z + l) 2 + p 2 + a|] 2mim 2 / z 2 + p 2 — 1 + za 2 (z — 1) 



a 2 i?3i?4 



(a 2 - 2i)i? 3 



+ c.c. 



7 



m 2 
4a 2 



+ 



2mim 2 (a 2 — 4) •m\ , m 2 / z 2 + p 2 — 1 + ia 2 (z — 1) 



(a| + 4)= 
m 2 [(z + l) 2 4-p 2 + a 2 
4a 2 i?3i?4 



(a 2 - 2*) 2 i? 3 



c.c. 



m\p 2 
2R\ 



(56) 



The metric is asymptotically flat which means that on the symmetry 
axis uo = 7 = when \z\ > 1. The balance of two Kerr particles then requires 
u) = and 7 = for \z\ < 1. The first of the latter two conditions yields, 
1 — z, i?3 = 1 + ia 2 + z, R4 = 1 — ia 2 + z, 

4/ 2 

2 [2mim 2 a 2 - z/(a 2 + 4)] = 0, (57) 
a 2 + 4 



after setting Ri 



whereas the second condition leads to 

4mim 2 / 2 (aJ — 4) 



0. 



(al + 4) 2 

From (JoTj) we get the value of v required for the equilibrium: 

2mim 2 a 2 



v 



a| + 4 
23 



(58) 



(59) 



On the other hand, equation (J5SJ) is satisfied when |a 2 | = 2. Introducing 
explicitly the separation parameter b, we finally obtain 

2b=\a 2 \. (60) 

Formula (|6U|) defines, for a given value of the angular momentum per unit 
mass of the superextreme Kerr constituent, the coordinate distance at which 
the equilibrium of two particles takes place. 

Therefore, this approximation scheme 3 too proves to be consistent with 
the exact double-Kerr solution. 

5. The Komar masses and angular momenta in the approximate 
solutions 

In this section we shall discuss the relation of the parameters rrii and of our 
approximation procedure to the genuine individual masses Mj and angular 
momenta Ji of the particles defined by the Komar integrals For the 

metric in the case of the cylindrical surface of integration we obtained 
the formulae [To] 




where Ui and U, i = 1,2, are the points of the z-axis which are taken as 
centers of the upper and lower bases of the cylinder enclosing the ith particle, 
Po being the radius of the bases. In our calculations we shall assume that the 
upper cylinder extends from z = to plus infinity (i.e., l\ — 0, U\ — +oo), 
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and the other one from minus infinity to z = (1% = — oo, u 2 = 0). The 
results of the evaluation of Mj and Jj can be summarized as follows. 

Scheme 1. Although in this scheme there are no equilibrium configura- 
tions and the particles are kept apart by a supporting strut, it is important 
for our analysis of the Komar quantities that the strut be massless in which 
case we can be sure that Mj and J, will not contain any contribution of the 
strut. Since the strut is massless when u = on it, we have to put v — 
in the formulae (J3*3"j) -(|37 JI . The integration with the aid of the formulae 1)61)) 
then yields 

Mi = Irrii, Ji = Pr/iidi. (62) 

Therefore, within the framework of this approximation scheme with v — 
the parameters m, represent the masses of the particles, and the parameters 
Oj coincide with the corresponding angular momenta per unit mass. 

Scheme 2. For the same reasons as in the previous case, in order to have 
correct estimations of Mj and Jj it is sufficient to consider a more general 
case of particles supported by a massless strut instead of considering config- 
urations of balancing particles. So, in the formulae ()44j ) -([46| ) we will take 
the constant v in the form (|4*9*|l to assure u = on the strut, and leave a,\ 
and a 2 as arbitrary parameters. The resulting expressions for mj and Jj are 

?Lf i 7/, ?2 4ajm 1 m 2 

Mj = /mj, Jj = Irriiai + / 7 rz -, (63) 

K + a 2 ) 2 + 4' v ; 

which means that in this scheme the parameters a, do not fully coincide with 
the Komar angular momenta per unit mass of the particles. 

Scheme 3. In the third approximation scheme the massless strut between 
the particles can be introduced by choosing v in the form ([59)1 . The evaluation 
of the Komar masses and angular momenta in (j61|) with the aid of (l55|) - (JoT)j) 
leads to the following result: 

2 2 4/77.17712(32 

Mj = /mj, J\ = I mxai, J2 — Iva^di + t — • (64) 

a 2 + 4 

Hence, the parameters mj in this scheme coincide with the Komar masses 
of the particles, and a\ with the angular momentum per unit mass of the 
subextreme constituent. The angular momentum per unit mass of the su- 
perextreme constituent is 

± = a 2 (l + l4^-), (65) 
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and it differs slightly from a 2 due to the presence of the second term in 
brackets. 

6. Discussion of the results 

Therefore, we have succeeded in working out correct approximate analogues 
of the exact double-Kerr solution by introducing a small parameter and 
expanding the Ernst complex potential in a series of powers of this param- 
eter. It turns out that the approximation procedure depends drastically on 
the qualitative relations between the parameters; it is not unique since a 
concrete physical situation dictates the form of the corresponding approxi- 
mation scheme. It is worthwhile mentioning that within all three schemes 
further variations of the procedure are possible. For instance, in scheme 1 it 
is possible to make the parameter v proportional to I 2 or I 3 . However, since 
these modifications cannot lead in principle to equilibrium configurations of 
two subextreme Kerr particles, we did not include them into this paper. At 
the same time, one should realize that by choosing v proportional to I 3 he 
commits himself to the calculation of the Ernst potential and correspond- 
ing metric functions up to the third order in I (the system of the balance 
conditions in this case takes the form 

2l 3 [m 1 m 2 {a 1 + a 2 ) - 2u] = 0, - mi m 2 l 2 = 0, (66) 

and implies vanishing of one of the masses.) 

A very important qualitative result with which the schemes 2 and 3 pro- 
vide us is the following: the location of a superextreme constituent in the 
approximate solution is not a point on the symmetry axis but it is a cut de- 
termined by two complex conjugate square roots, say, Ri = R 2 . One should 
be aware of this non-evident result when he wants to develop an approxima- 
tion procedure for spinning particles without relating it to a concrete exact 
solution. 

In view of the above said and of our Section 2 it would be of interest to 
compare Bonnor's solution (|1J with our approximate solution from scheme 1. 
The approximation scheme employed by Bonnor in ^3] treats the parame- 
ters rrii and ji as small parameters of the same magnitude in I. This means 
in particular that the quantities a« = ji/ m * m his approach are not small 
constants and, consequently, both spinning particles are superextreme be- 
cause a 2 > m 2 . The fact that the superextreme constituents of the solution 
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are placed at the points ±6 of the symmetry if they were subex- 

treme constituents makes Bonnor's approximation procedure contradictory 
from the point of view of the double-Kerr solution. For the solution (jjj) 
to represent two subextreme Kerr particles it is necessary to introduce the 
small parameters aj into (@J) by additionally changing ji to Im^ai (to make 
ji proportional to I 2 ) and eliminate all the arising terms of the third order 
in I. Then after such an extensive make-up the resulting solution will be a 
particular specialization [y = 0) of our approximate solution from scheme 1 
(the lack of one arbitrary parameter in the solution Q was discussed in [28J). 
Mention that in the static limit the solutions (HJ) and (|35|l -(|37 j) coincide and 
represent correctly an approximation to the double-Schwarschild spacetime. 

Turning back to our approximate solutions from schemes 2 and 3, below 
we would like to demonstrate that the condition (|5Uj) determining the balance 
of two superextreme Kerr particles and the condition (j6T?j) determining the 
balance of a sub- and a superextreme Kerr constituents agree with the exact 
results on the double-Kerr equilibrium problem not only qualitatively but 
also quantitatively Indeed, in our paper |TH! we obtained the general relation 
for balancing Kerr particles 



where J = M X A X + M 2 A 2 , M = M 1 + M 2 , M; are Komar masses, J { = MiAi 
are Komar angular momenta, and s is the coordinate distance separating the 
particles (s = 2b in the zeroth order in I). 

First let us assume that both constituents in (|67|) are superextreme. Then, 
substituting Ml — > IM\, M 2 — > IM 2 into (|57j) and restricting ourselves to only 
zeroth order in I we obtain 



This is exactly the condition (|5Uj). 

Assuming now that the first particle is subextreme and the second one is 
superextreme, we substitute Mi — > I Mi, M 2 — > IM 2 , A\ — > IA\ into (joT|) and 
get, after retaining only the zero-order terms in I: 



J±(M + s) 2 + s(A 1 + A 2 )=0, 



(67) 



± s 2 + s(A 1 + A 2 ) = s = \A 1 + A 2 



(68) 



± s 2 + sA 2 = s=\A 



2 h 



(69) 



which is precisely the condition ([BOJl. 
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7. Conclusions 

The formalism developed in this paper opens possibilities for the unified 
treatment of the exact and approximate solutions of Einstein's equations. It 
is essentially based on Sibgatullin's approach to the construction of exact 
solutions from the prescribed axis data that permits establishing very strict 
and clear connections between the parameters of the solution and the corre- 
sponding physical properties they determine uniquely through the relativistic 
multipole moments. 

We have constructed approximate analogues of the well-known exact 
double-Kerr solution and demonstrated that they permit equilibrium con- 
figurations involving superextreme particles. In the case of the double-Kerr 
equilibrium problem the exact approach looks by far more attractive since 
the balance conditions can be solved analytically in a unified form, whereas 
a specific approximation procedure has to be worked out for each concrete 
type of the Kerr particles. The situation, however, may change in the case 
of more complex exact solutions describing the many-body configurations 
for which the balance conditions cannot be solved analytically in general. 
Then it is anticipated that the approximate procedure discussed in this pa- 
per may become the most efficient tool for obtaining the general relations for 
balancing particles, especially in the presence of the electromagnetic fields. 
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Figure 1: Exact analogue of the approximate solution (4): a specific four- 
body system formed by two compound objects. 
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Figure 2: Three principle configurations of two Kerr particles treated by 
approximation schemes 1-3. 
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